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Abstract 

We use AdS/CFT duality to study the thermodynamics of a strongly coupled J\f = 2 super- 
symmetric large N c SU(N C ) gauge theory with Nf = 2 fundamental hypermultiplets. At finite 
temperature T and isospin chemical potential /i, a potential on the Higgs branch is generated, 
corresponding to a potential on the moduli space of instantons in the AdS description. For 
/i = 0, there is a known first order phase transition around a critical temperature T c . We find 
that the Higgs VEV is a suitable order parameter for this transition; for T > T c , the theory is 
driven to a non-trivial point on the Higgs branch. For /i ^ and T = 0, the Higgs potential 
is unbounded from below, leading to an instability of the field theory due to Bose-Einstein 
condensation. 
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1 Introduction 



In its original form, AdS/CFT duality 01210] relates theories of closed strings in asymptotically 
AdS spaces to large N c gauge theories with matter in the adjoint representation. Fields in 
the fundamental representation may be added by including an open string sector through the 
introduction of branes probing the supergravity background. Much effort has gone into studying 
dualities of this type, motivated largely by the goal of finding a supergravity background dual 
to QCD. 

The first example of AdS / CFT duality for a theory with fundamental representations related 
a conformal M = 2 Sp(N) gauge theory to string theory in AdS^ x S 5 /Z 2 , with D7-branes 
wrapping the Z 2 fixed surface with geometry AdS§ x S 3 |H E] • A different approach to flavour 
in AdS/CFT has been considered in 0. In jjj, the duality of |HE] was extended to an J\f = 2 
SU(N C ) theory with Nj massive fundamental hypermultiplets, essentially by removing the Z2 
orientifold, which was justified by the fact that a probe D7-brane wrapping a contractible S 3 does 
not lead to a tadpole requiring cancellation. Although the field theory is not asymptotically free, 
it has a UV fixed point in the strict N c — > 00 limit. Following this, there have been a number of 
papers generalizing the duality to confining theories with fundamental representations [S] -|19j. 
including non-supersymmetric examples in which spontaneous chiral symmetry breaking by a 
ipip quark condensate occurs [20] -[2H]- 

At zero temperature and vanishing quark mass, the M = 2 gauge theory has a non-trivial 
Higgs branch, i.e. a moduli space of vacua on which the scalar components of the fundamental 
hypermultiplets have expectation values. For finite quark mass, the moduli space includes a 
mixed Coulomb-Higgs branch: the fundamental hypermultiplets may have non-zero expectation 
values if the vector multiplet scalars have particular VEV's equal to the quark mass. For both 
the massless and massive case, the moduli space of field theory is described in the AdS picture 
by instanton configurations on the D7-branes (HHl Ell lEH lEHj- Self dual field strengths are 
solutions of the D7-brane equations of motion due to a conspiracy between the Yang-Mills 
and Wess-Zumino terms in the D7-brane action. In general, this conspiracy is destroyed by 
supersymmetry breaking deformations, which give rise to a potential on the moduli space of 
instantons, corresponding to an effective potential on the Higgs branch. We will compute 
this potential at finite isospin chemical potential and temperature. We focus on the potential 
generated on a slice of the moduli space corresponding to a single instanton centered at the 
origin. This slice is parameterized by the instanton size, which is dual to a particular Higgs 
VEV. 

At finite isospin chemical potential, we find the expected result that the theory is destabilized 
via an effective negative mass squared for the moduli, leading to Bose-Einstein condensation. 
The negative mass squared term in the effective potential is related to the metric on the Higgs 
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branch, which is correctly reproduced by the dynamics of spinning instantons on the D7-branes. 
To stabilize the theory, a positive scalar mass which is larger than the chemical potential would 
need to be introduced. 

Next we consider finite temperature deformations. Some finite temperature properties of 
theories with fundamental representations have been studied previously using AdS / CFT duality 
PIH EH]- As shown in j20J there is a finite temperature first order phase transition as 
the ratio of the temperature to the quark mass is varied in the M = 2 theory of [7] (a similar 
transition was discussed in |22j)- In the AdS description, this transition corresponds to a change 
in the topology of the D7-brane embedding in an AdS-Schwarzschild background. 

We explore the finite temperature behavior of this theory in more detail by computing the 
effective potential generated on the Higgs component of the moduli space. We find that the 
Higgs expectation value is an order parameter for the first order phase transition. At < T < T c 
and fi = 0, the Higgs VEV is driven to the origin of moduli space. However for T > T c we find 
a surprise; the instanton size is driven towards a non-zero value, suggesting the existence of a 
vacuum in which the theory is higgsed. This is in contrast to the weak coupling behavior of the 
theory, for which the one-loop finite temperature effective potential implies that the origin of 
moduli space is at least metastable. 

The organisation of this paper is as follows. In section 2, we review the AdS description 
[3*U| ITTTj of the mixed Coulomb-Higgs branch in the Karch-Katz M = 2 theory In section 
3, we consider the same theory at finite chemical potential. Section 4 is devolved to the finite 
temperature case. In section 5 we briefly summarize our results and discuss future developments. 

2 SUGRA dual of an J\f = 2 theory with fundamental 
representations 

We consider an N = 2 gauge theory which is dual j7j to string theory in AdS$ x S 5 with Nf D7- 
branes wrapping a surface which is asymptotically AdS§ x S 3 . The matter content of this gauge 
theory is that of the M = 4 SU(N C ) gauge theory together with Nf massive hypermultiplets in 
the fundamental representation. In M = 1 superspace, the Lagrangian is 




(i) 



where the superpotential W is 



(2) 
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Coordinates 


1 2 3 4 5 6 7 


8 9 


D3 




D7 




x^ y m (or y, Q 3 ) 







Table 1: Index conventions 



The superfields Q l and Qi, labeled by the flavor index % = 1---Nf, comprise the M = 2 
fundamental hypermultiplets. 

This theory is not asymptotically free and, at finite N c , the corresponding string background 
suffers from an uncancelled tadpole. However, as in [Zj, we focus strictly on the N c — > oo limit 
with fixed Nf. In this case there is a non-trivial UV fixed point for the 't Hooft coupling, 
while the dual AdS string background does not suffer from a tadpole problem since the probe 
D7-branes wrap a contractible S 3 . 

In coordinates which will be convenient for our purposes, the AdS$ x S 5 background is 

r 2 R 2 

d s 2 = —dx»dx» + — (dy 2 + y 2 dn 2 3 + dZ l dZ l ) , (3) 

i=l 

$ 

e = 9s, 

4 

F (5) = dC (4 ) = 4i? 4 (F s5 +* V s s), C (4 ) 1 0123 = ^ c?x° A dx 1 A c/x 2 A rfx 3 , 
R A = Ang s N c a' 2 , 

where V55 is the volume form on S* 5 , and r 2 = y 2 + Z t Z t . The 't Hooft coupling in the dual 
gauge theory is A = g 2 N c = g s N c . 

The fundamental hypermultiplets arise from Nf D7-branes embedded in this geometry. Their 
action is 

Sm = ~T 7 J d t xd*yy/-detG ab \ PB = - T 7 j d*x dy dtt 3 y 3 \/ 1 + (Z^) 2 , (4) 

where Gpb is the pull-back of the metric, a prime indicates a derivative with respect to y, and 
we have assumed an embedding independent of the coordinates x as well as the coordinates on 
S 3 . There is a manifest minimum of the action when Z[ — 0. We will choose 

Z 1 = m, Z 2 = 0, (5) 
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other solutions being related by the U(1)r symmetry corresponding to rotations in the plane 
spanned by the Z l . The induced metric is 



R 2 



{dy 2 + y 2 dQ 



3)1 



(6) 



with r 2 = y 2 + m 2 . The parameter m corresponds to the mass of the fundamental hypermul- 
tiplets. For m — 0, the geometry is AdS§ x S 3 , while for m ^ 0, the geometry approaches 
AdS 5 x S 3 at large r. The 5* 3 component of the D7-geometry contracts to zero size at r = m. 



So long as Nf is held fixed in the limit N c — > 00 with fixed A 
back-reaction of the D7-branes on the bulk geometry. 



g s iV c 3> 1, one can neglect the 



2.1 The Higgs branch 

When the theory has massless quarks, the fundamental scalars q l and (denoting bottom 
components of chiral superfields by lowercase letters,) have non-zero expectation values on 
the Higgs branch 1 while the adjoint scalar 03 of the TV = 2 vector multiplet vanishes. For 
non-zero mass, m, and vanishing 3 , the fundamental hypermultiplets are massive and a pure 
Higgs branch does not exist. However there is a mixed Coulomb-Higgs branch when 3 has an 
expectation value such that some components of the hypermultiplets are massless. An example 
of a point on a mixed Coulomb-Higgs branch is given by a diagonal 03 for which all but the last 
k entries are vanishing, 



/o 



\ 



-m 



(7) 



-m) 



lr The scalars tfii and 02 belonging to the adjoint hypermultiplet may also have non-zero expectation values 
on the Higgs branch. 
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In this case, the F-flatness equations qi{<pz + m ) — (03 + rn)q l = permit fundamental hyper- 
multiplet expectation values in which only the last k entries of q l and q~i are non-zero, 



(0 



fi- 



ll 



Pki) 



(8) 



w 



There are additional F and D-flatness constraints which we have not explicitly written. 

In string theory, nonzero entries in (JHJ) physically correspond to D3-branes which are coinci- 
dent with and dissolved within the D7-branes. Dissolved D3-branes can be viewed as instantons 
in the eight- dimensional world- volume theory on the D7-branes JHE], due to the Wess-Zumino 
coupling 

S wz = ^(2W) 2 / C\% A tr(F A F) . (9) 
9s J 

There is a one to one map between the moduli space of Yang-Mills instantons and the Higgs 
branch of this M = 2 theory The ADHM constraints from which instantons are constructed 
[3*7j are equivalent to the F and D-flatness equations JHE1IHH] (see also J20J for a review). 



2.2 Supergravity description of the Higgs branch 

Because of the known one-one correspondence between instantons and the Higgs branch, one 
expects that instantons solve the equations of motion of the non-Abelian 2 Dirac-Born-Infeld 
action describing D7-branes embedded in (JHJ) according to (JHJ). The existence such solutions is 
a non-trivial consequence of AdS / CFT duality JHD1 EH HI] • 

The effective action describing D7-branes in the AdS background (JHJ) is 



S = - T 7 (2W) 2 ( -— J rf 8 eC (4) A tr(F A F) + J d 8 £ — \/ - det G Tr (F a(3 F a ^) ) + 



(10) 



where we have not written terms involving fermions and scalars. This action is the sum of a 
Wess-Zumino term, a Yang-Mills term, and an infinite number of corrections at higher orders 
in a' indicated by • • • in IjlOj) . The correspondence between instantons and the Higgs branch 
suggests that the equations of motion should be solved by field strengths which are self-dual 
with respect to a flat four- dimensional metric. 



2 The existence of a Higgs branch requires at least two flavors, or two D7-branes. 



6 



In this paper, we work to leading order only in the large 't Hooft coupling expansion gener- 
ated by AdS / CFT duality, which allows us to neglect the higher order terms in the a' expansion 3 
of the action. Constraints on unknown higher order terms arising from the existence of instanton 
solutions, as well as the exactly known metric on the Higgs branch, were discussed in pltH KT] . 

The induced metric (JBJ) can be written as 



ds 2 = ^dx»dx» + ^rY dy m dy m , (11] 
R z r z ^— ' 

m=l 

with r 2 = y m y m +m 2 . Field strengths which are self dual with respect to the flat four-dimensional 
metric ds 2 = Ylm=i dy m dy m solve the equations of motion, due to a conspiracy between the 
Wess-Zumino and Yang-Mills term. Inserting the explicit AdS background values (jHJ) for the 
metric and Ramond-Ramond four- form, into the action (jl(J|) for D7-branes embedded according 
to ((SJ), with non-trivial field strengths only in the y m directions, gives 

/ r r 4 / i \ 

S TjylTtOL ) / d xd y — . I "Z^-mnrsF m nF rs ~\~ F mn F mn ) 

J Ag s K \ 2 J ^ 

= -T 7 (A7ra') 2 [ d 4 xd 4 y- 



where F~ n = \(F mn — \e mnrs F rs ) . Field strengths F~ n = 0, which are self-dual with respect to 
the flat metric dy m dy m , manifestly solve the equations of motion 4 . These solutions correspond 
to points on the Higgs branch of the dual M = 2 theory. If m ^ this is, strictly speaking, 
a point on the mixed Coulomb-Higgs branch, with expectation values of the form (jJJ), (|5|) . 
We emphasize that in order to neglect the back-reaction due to dissolved D3-branes, we are 
considering a portion of the moduli space for which the instanton number k is fixed in the large 
N c limit. 



2.3 A slice of the Higgs branch 

For simplicity, we consider the case Nf = 2, which is the minimum value for which a non-trivial 
Higgs branch exists. We will focus on a slice of the Higgs branch (or mixed Coulomb-Higgs 
branch) corresponding to a single instanton centered at the origin, y m = 0. 

In "singular gauge", the SU(2) instanton (using the same conventions as jlOj) is given by 



2Q O nmHn 

y 2 (y 2 + Q 2 Y 



3 The dimensionlcss expansion parameter is a' / R 2 = l/y/X. 

4 Anti-instantons, with F + = 0, correspond to non-supersymmetric configurations which do not solve the 
equations of motion. 
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where Q is the instanton size, and 



On 



O n O r , 



0„ 



. \OmOri <T n (J r , 



O r , 



(if A 



2x2; 



Or, 



O 



In = (-ir,l 2X 2) 



(14) 



with t being the three Pauli-matrices. We choose singular gauge, as opposed to the regular 
gauge in which A n = 2a mn y m / (y 2 + Q 2 ), because of the improved asymptotic behavior at large 
y. In the AdS setting, the Higgs branch should correspond to a normalizable deformation of 
the AdS background at the origin of the moduli space. We see that the solution falls to zero 
at large y and the deformation in the interior is controlled by the dimension two parameter Q 2 
which corresponds to the VEV of q l qi. The choice of singular gauge will be of particular use 
when computing the effects of the chemical potential. The singularity of (|13|) at y m = is not 
problematic for computations of physical (gauge invariant) quantities. 

For m = and zero instanton number, the background geometry preserves the symmetries 
50(2,4) x SU(2) L x SU(2) R x U(1) r x SU(2) f . The 50(2,4) isometry of the AdS 5 factor 
corresponds to the conformal symmetry of the dual gauge theory. The SU(2)l x SU(2)r ~ 
50(4) acts in the obvious way on the coordinates y m . The SU(2)l factor corresponds to a global 
symmetry, while SU(2)r and U(1)r, which acts on (Z 1 ^ 2 ), correspond to R-symmetries 5 . 
SU(2)f is a gauge transformation on the D7-branes which is constant at the boundary of AdS. 
The action at the boundary of AdS corresponds to the flavor symmetry of the dual gauge theory. 

In the presence of the instanton ([13)1 . and for m^O, the symmetries are broken to 50(1, 3) x 
SU (2) l x diag(5?7 {2)r x SU (2) /) corresponding to a point on the mixed Coulomb-Higgs branch 

/ o \ /o \ /o \ 





\tiaQ J 



\ 







V 



-m 



(15) 



/ 



where qi a are the scalar components of the fundamental hypermultiplets, with flavor index 
i = 1, 2 and SU(2)r index a = 1, 2. The scalars <p a belong to the adjoint hypermultiplet, while 
03 is the adjoint scalar belonging to the vector- multiplet. 

In general, one expects that supersymmetry breaking deformations, such as a finite chemical 
potential or temperature, lift the vacuum degeneracy of the Higgs and Coulomb branch. In 
the dual supergravity description, the lifting of the Higgs branch at finite temperature occurs 
because the conspiracy between the Wess-Zumino and Yang-Mills terms, necessary for self- 
dual field strengths to be solutions, no longer occurs. The potential which is generated on the 
Higgs branch can be computed by evaluating the D7-brane action on the space of self-dual field 



5 There is no apparent distinction between SU(2)l and SU(2)r in the geometry of the AdS background. The 
distinction arises from the gauge field couplings to the RR four-form. 
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strengths. In the following we will compute the potential V(Q) generated on the slice of the 
Higgs branch dual to the single instanton 

3 Finite chemical potential, zero temperature 

As the simplest example of a potential generated on the Higgs branch, we first consider the case 
of finite chemical potential and zero temperature. For two flavors, there is a U(2) flavor symme- 
try, with Cartan generators 1 and 03, corresponding to baryon number and isospin respectively. 
We will consider a nonzero chemical potential for the isospin 6 . To include the chemical potential 
|%2~j we allow a spurious gauge field associated with the r 3 component of isospin to acquire a 
VEV, fi, in its A component. This includes generic fermion and scalar Lagrangian terms for 
fields with isospin charge e of the form 

5C = -/ie^r 3 7 V + /i 2 e 2 |0| 2 . (16) 

The first term is a source for the fermionic isospin number density. In the path integral, this 
term places the theory at finite density. The second term is an unbounded scalar potential 
which renders the theory unstable, such that Bose-Einstein condensation is expected. We will 
reproduce this run-away behaviour of the scalar potential using the AdS/CFT description. 
Equivalently, we may set A T = by a gauge transformation e lfieT , and perform the path integral 
with boundary conditions corresponding to a spinning configuration space. 

AdS/CFT duality relates global symmetries of the boundary theory to gauge symmetries 
in the bulk supergravity. The prescription for turning on a chemical potential in the AdS 
description is to turn on a background (non-normalizable) flat gauge connection A T = fi for 
the associated gauge symmetry. Equivalently, one may consider a spinning AdS background. 
Examples of the AdS description of a chemical potential in various contexts have appeared in 

E3-E3- 

The effect of the chemical potential on the Higgs branch can be studied in the AdS description 
by computing the action for rigidly rotating instantons with moduli A4 l , 

A n = e ^3 ^instanton ^m^ ^iy-i^z _ ^ 

Equivalently, we may add a background A in the fixed instanton background, 

A = , A n = ^ nstanton . (18) 

On the slice of the Higgs branch corresponding to the single instanton configurations (|13|) with 
modulus Q, the effective potential at quadratic order in \i can be determined by inserting ()18|) 

6 Turning on a chemical potential for the baryon number has no apparent effect here since none of the fields 
on the D7 brane world volume are charged under this symmetry. This was pointed out to us by A. Karch. 
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into the D7-brane action J d 4 x V(Q) = —Sdi giving 

V(Q) = T7 { ^^ J d 4 ytr (\ iy2 + R f )2 F- n F~ n + 2F mM FwT (19) 

+ (y 2 + m 2 ) 2 ^ a(3V V J ' 

with y 2 = y m y m . We have split the action into the pieces involving F in the x and y directions, 
indicated by Greek and Roman indices respectively, as well as mixed terms. For the background 
(|18|) . the only non-zero contribution to the potential comes from the mixed term tr F /im F um rj fJ ' u = 
-tr [A ,A n ] 2 , giving 

v(n\ t 2 ( W ) 2 2 f m Q* rp 2(4ttV) 2 
V(Q) = -T 7 /i / d y = -T 7 p Q . (20) 



9s J y 2 (y 2 + Q 2 ) 2 9s 

Note that the term ^fgF^ m F^ m in the D7-action is also the term which determines the 
metric (two-derivative effective action) on the Higgs branch. To quadratic order in the chemical 
potential, the result pUj) is exact. Higher dimension operators in the DBI action with two Greek 
indices, which have been omitted from (|19|). must vanish in the instanton background, otherwise 
there would be corrections to the Higgs branch metric in a strong 't Hooft coupling expansion 
[3?H I3T] . By the same token, the term in the potential on the Higgs branch which is quadratic 
in the chemical potential is given exactly by the F iim F lim term. 

In light of (|20p. the M = 2 theory is unstable at finite chemical potential. Physically, the 
instability is due to Bose-Einstein condensation. This instability will not be cured by taking 
into account the higher derivative terms in the effective action, whose contribution does not 
effect the large Q asymptotics. Potentially, the instability could be removed by the inclusion of 
an appropriately large mass term for the scalars, either as an explicit additional perturbation, 
or generated radiatively from the inclusion of some other relevant perturbation. The finite 
temperature configurations studied below, however, will not be able to cure this instability 
since - as shown below - the relevant deformation, T 4 , has dimension four and can at best 
generate a potential T 8 /Q A . This does not stabilize the theory. 

4 Finite temperature, zero chemical potential 

Before investigating the large 't Hooft coupling behavior of the effective potential on the Higgs 
branch, let us review the situation to one loop in perturbation theory. In a supersymmetric 
theory, the one-loop effective potential is 

d 2 W 



hi 



(21) 
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where V is the zero temperature effective potential, W is the superpotential, and 0j are the 
scalar components of chiral superfields. For the theory we are considering, this potential drives 
the system towards the origin of moduli space. In fact it is often said that supersymmetric 
theories with flat directions do not have symmetry breaking vacua at finite temperature jlH] 
(for some exceptions, see [SniO])- However, for large 't Hooft coupling, we will find that the 
theory we consider does indeed exhibit symmetry breaking at high temperature. The origin of 
the Higgs branch becomes a maximum of the effective potential for T > T c . 

The gravitational dual of M = 4 gauge theory at large 't Hooft coupling and finite tem- 
perature is obtained by replacing the background (JHJ) with the AdS-Schwarzschild black-hole 
background The latter belongs to a general class of supergravity solutions which, in a 

choice of coordinates convenient for our purposes, have the form 

4 2 

ds 2 = f(r)(dx 2 + g(r)dr 2 ) + h(r)(%2 dy m dy m + y ^dZ i dZ i ), 



e 



m=l i=l 



F (5) = ^(Vss +* V s b) = dC {4) , C (4 )|oi23 = s(r) dx° A dx 1 A dx 2 A dx 3 , 

r 2 = y m y m + Z { Z\ (22) 
For the AdS-Schwarzschild solution, we have 

f(r) = W , »W = (j^ ) , Mr) = ^, ^r) = - (l + — j . (23) 

The coordinates x are the spatial coordinates of the dual gauge theory and r is the Euclidean 
time direction, which is compactified on a circle of radius b' 1 , corresponding to the inverse 
temperature. Note that the temperature T ~ b only enters to the fourth power. 

The D7 embedding in this background to describe our TV = 2 gauge theory with quarks is 

Z 2 = 0, Z x = z{y), (24) 

where y 2 = y m y m . The induced metric G| p b takes the form 

ds 2 m = f(r) {dx 2 + g{r)dr 2 ) + h(r) (dy m dy m + ' (y m dy m A , (25) 

with r 2 = y 2 + z 2 (y). To lowest order in a', the D7-brane action is 

-T 7 J d 8 £ e"*v /:r detG| pb = -2tt 2 T 7 J d A x dy y s ''<f>(r) f(r) 2 \/gJ^ h(r) 2 Jl + z'{y)\ (26) 

and the specific embeddings z(y) can be found by solving the Euler-Lagrange equation with 
boundary conditions 2(00) = m, z'(oo) = 0. 
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At large y the geometry returns to AdS where the asymptotic embedding solution (J24j) takes 
the form z(y) ~ m + c/y 2 . Given m, the parameter c, corresponding to a tpip condensate, 
is determined by requiring the solution to be smooth and normalizable. These solutions were 
studied in detail in [2U] . 

The potential generated on the Higgs branch is obtained by calculating a contribution to the 
D7-brane action involving the gauge field F a ^ at second order in a'. Specifically, one evaluates 
the action on the space of field strengths which are self-dual 7 with respect to the induced metric 
in the directions transverse to r, x\ 



(27) 



V = T7(2 ^ )2 (±- J d A y C$ 23 e mnrs Tr F mn F rs -I J d A y e~* v^detG Tr F mn F mn ^j 

where F mn is self-dual with respect to the metric 

ds\ = h{r) ((1 + z\y) 2 )dy 2 + y 2 dQj) . (28) 

This metric is conformally flat. With new coordinates y(y) such that ds 2 = a(y)(dy 2 + y 2 dQ 2 ), 
the instanton configurations (self-dual field strengths) take the usual form. 

In the 6 = (zero temperature) limit, the analysis of section (|2.2j) holds and the potential 
(J27|) becomes flat. At finite temperature, the conspiracy between the Wess-Zumino term and the 
Yang-Mills terms fails, giving a potential on the moduli space of instantons. We now compute 
the form of the potential for the slice of moduli space corresponding to an instanton centered 
at the origin. 

4.1 Computation of the Higgs Effective Potential 

We wish to evaluate the D7-brane action on the space of fields strengths which are self-dual 
with respect to the transverse part of the metric 

ds{ = h{r) I dy m dy m + ( ^M) (y m dy m ^ 2 



y 

This metric is conformally flat. There exists a coordinate transformation 

y m = Ay)y m (29) 

such that the induced metric has the form 

ds 2 D7 = f(r) {dx 2 + g{r)dr 2 ) + ^dy m dy m . (30) 



7 Thcrc are couplings between world-volume scalars and field strengths at higher orders in a' which could 
alter the embedding. However we only consider the leading term in a large 't Hooft coupling expansion for which 
these couplings can be neglected. 
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The function J(y) is obtained by solving 



J'V 2 J 1 1 

7 







(31) 



y J y 2 

subject to the boundary condition J(y — > oo) = 1. In the coordinates y m , the single instanton 
centered at the origin has the usual form, 

2Q 2 a nm y n 



Am y 2 (y 2 + Q 2 )' 

The embedding function z(y) is determined by shooting techniques as in 
the Higgs branch is then given by Sd7 = — J d A xV(Q), or 

(27ra') 2 



(32) 

The potential on 



V(Q) = T 7 



2 
4& 



1 y d 4 yC$ 23 e mnrs TrF mn F rs -~ J d 4 y e^V^tGTrF mn F mn ^j 
d 4 y (s(r) - 0(r)/(r) 2 v^(r) 



96 g 4 
(y 2 + g 2 ) 4 



(33) 



70/(5-3) J y 3 dy (s(r)-0(r)/(r) 2 V^ 



96 g 4 j(i/) 4 
(J(y)V + Q 2 ) 4 



</(y) 



For the AdS-Schwarzschild background 
6(2ttV) 2 



V{Q) 



T 7 



g s R 4 



dy 



we obtain 
bW Q 4 J{y) 4 



( J(y) V + %) 2 ) 2 ( J(y) v + g 



J 



(34) 



4.2 Large VEV Asymptotics 



Without resorting to numerics, some qualitative features of V{Q) at large Q (or large instanton 
size) can be determined from the large y behavior of the induced metric and RR-form given in 
1 together with the asymptotic behaviour of the embedding, 



z(y) = m + — + 



y- 



(35) 



At large y we can neglect the condensate and the leading terms in the potential will depend 
on m and on the deformation parameter of the background b. The D7-brane action for a large 
instanton, which has support at large y is 



d 4 xV{Q) = -S D7 



T-, 



6(2vr 2 a 
9sR A 

6(2ttV) 2 

g s R A 



d 4 x / dy y 3 



(y 2 + Q 2 ) 4 [y 



b s 2b s m 2 



If 



+ ... 



(36) 
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This asymptotic approximation is valid when y > Q but gives IR divergent integrals. We have 
evaluated the integrals down to y/Q — 1. In the next section we evaluate the full potential 
including the correct, convergent, IR behaviour. 

The asymptotic potential terms are simply what one would expect on dimensional grounds. 
The potential must vanish as b — > and the geometry becomes AdS. Since b enters as b A , the 
first Q dependent term must take the form b 8 /Q 4 . 

4.3 Numerical computation of V(Q) 

To compute V(Q) in general requires knowledge of the embedding function z(y), which has 
been computed by a numerical shooting technique in [20]. Imposing boundary conditions cor- 
responding to the large y behaviour (|35|). and requiring smooth behaviour in the interior, such 
that an RG flow interpretation is possible, leads to a dependence c(m, b) of the chiral quark 
condensate (ipip) on the quark mass and on the temperature. Depending on the ratio m/b, 
there are two types of solutions, which differ by the topology of the D7-branes. At large r 
(or y) the geometry of the D7-branes is AdS§ x S 3 and the topology of the r — ► oo boundary 
is S 1 x R 3 x S 3 . For sufficiently large m/b, the S 3 component of the D7-geometry contracts 
to zero size at finite r > b. In this case the D7-brane "ends" before reaching the horizon at 
r = b. However, for sufficiently small m/b, the D7-brane ends at the horizon, at which point the 
thermal S* 1 contracts to zero size. Both these types of solutions are plotted in figure ^ There is 
a first order phase transition at the critical value of m/b ~ 0.92 where the two types of solution 
meet |2*Ul I35j . The ("tpip) condensate is non-zero on both sides of this transition, although there 
is a discontinuous jump in its value. 



z(y) 




m=0.92 
m=0.91 



Figure 1: Brane embeddings in the AdS Schwarzschild background for different values of the 
quark mass. In the plot we have set b = 1. 
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Using the machinery discussed above, we can study the behaviour of the scalar VEV 
parametrizing the Higgs branch on either side of the transition. We evaluate the potential 
()27I34|) numerically, with a single instanton centred at the origin for different values of the ratio 
m/b. The results are plotted in figure 121 




V(Q) V(Q) 



Figure 2: Potential V(Q) as a function of the instanton size / Higgs VEV Q for various values 
of the quark mass m (we set b = 1 here); in fig. a) we sample the whole spectrum from m = 
to m — > oo, in fig. b) we vary m close to the phase transition region. The potential for m — > oo 
is flat and coincides with the horizontal axis. 



Figure 121 shows that there are two phases for the scalar Higgs VEV matching the two phases 
(topologies) of the D7-brane embedding. These phases are divided by the first order phase 
transition at a critical value of m/b, where the parameter b of the AdS/Schwarzchild background 
is proportional to the temperature. In figure El we have chosen units such that b = 1. For 
m/b — > oo, the D7 branes probe the AdS-like part of the geometry, such that the potential is 
flat (dotted line on the left hand side of figure |2J). For all other values of m/b, the potential 
approaches this constant asymptotic value from below for large Q. When the value of m/b is 
decreased, starting at m/b — > oo, a minimum of the potential forms at Q — 0, with increasing 
depth. At the critical value m c /b ~ 0.92, a new solution for the brane embedding appears, for 
which the potential has a minimum at Q ^ 0. For m/b — > 0, V has a maximum at Q = and 
a minimum of decreasing depth for Q ^ 0. At m/b = (solid line in the plot) there is still a 
shallow minimum of V at Q ~ 0.856. Thus, the Higgs VEV is an order parameter for the first 
order phase transition. We emphasize that the symmetry breaking phase occurs at small m/b, 
or high temperature, contrary to the usual expectation. We plot the value of the Higgs VEV as 
a function of m/b in Figure 01 
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0.9125 0.915 0.9175 0.92 0.9225 0.925 



m q /T 



Figure 3: Position of the minimum of the potential Qq versus the bare quark mass m q , zoom of 
the critical region. 

Note that the expectation value Qq is a multivalued function oim/b in a region near the first 
order critical point. There are several regular embeddings in this region, with the true vacuum 
corresponding to the solution with lower free energy. 

Thus we find that for T > T c , the theory is driven to a Higgs phase or, more precisely, a point 
on the mixed Coulomb-Higgs branch. We have not determined whether or not this vacuum is 
merely metastable, which would be the case if the origin of the Coulomb branch has lower free 
energy. 

Since the less supersymmetric phase occurs for T > T c , the phase transition discussed here is 
quite different from the deconfinement transition. Nevertheless it is interesting to note that our 
phase transition is first order, as is found in lattice simulations of the deconfinement transition 
in large N SU(N) gauge theories [52] . 

4.4 Geometric aspects of the transition 

It is worth stressing the link between the transition of the Higgs VEV between the two phases 
and the change in the embedding topology of the D7 brane. From figures 1 and 2, we see that 
the Higgs potential only has a non-zero minimum when the D7-branes reach the horizon. There 
is a repulsive effect in the vicinity of the horizon which causes the instanton to expand. In 
figure E] we plot the value of the Higgs VEV versus the minimum value of the coordinate y to 
which the D7 branes extend. When the D7-branes do not reach the horizon, the S 3 component 
of their geometry contracts to zero size at the endpoint, which occurs at y — 0. On the other 
hand, when the D7-branes reach the horizon, the S 3 component of their geometry remains finite 
size and y ^ at the endpoint. The dependence of the VEV on the minimum value of y is 
almost a straight line with slope roughly equal to 1.25. The solutions that reach the horizon are 
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distributed along this line, while the ones which end before reaching the horizon accumulate at 
the origin. 

The phase transition region (m/b ~ 0.91) has some additional interesting structure. In this 
region there are three regular solutions for the D7 embedding for each value of m q (see 122])- 
Some of them end on the horizon and some do not. The three flows are shown for a particular 
case in fig. Ek>. These different embeddings give different predictions for the Higgs VEV. We also 
show the potential for Q on each of these embeddings in figure and the minimum vs m can be 
found in figure Eb- There are three values of Q supplied in this transition region. Clearly there 
will be a first order transition where the embedding jumps from solutions ending on and off the 
horizon as m is changed, corresponding to a discontinuous transition in the value of the scalar 
VEV. 

5 Conclusions 

We have investigated the thermodynamics of a large N c strongly coupled M = 2 gauge theory, 
using AdS/CFT duality to compute the potential which is generated on the Higgs branch at 
finite temperature and isospin chemical potential. In the AdS description, this potential is a 
potential on the moduli space of instantons. 

We have shown non-perturbatively that the chemical potential destabilizes the Higgs poten- 
tial, giving rise to Bose-Einstein condensation. The 0(fi 2 ) term in the potential can be computed 
exactly, as a consequence of the non-renormalization of the metric on the Higgs branch, which 
is correctly reproduced by instanton dynamics in the AdS description. The instability due to 
the chemical potential is not cured at finite temperature. 

It would be very interesting to introduce a chemical potential for Baryon number in the AdS 
setting. Unlike the chemical potential for isospin, this is known to be very difficult to study 
on the lattice due to a complex fermionic determinant. In the AdS setting, it is also easier to 
introduce a chemical potential for the isospin, as we have done here, although the problems 
arising for baryon number are probably much more tractable than those which arise on the 
lattice. 
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Figure 4: Position of the endpoint y m in of the embedded brane on the horizon versus the 
potential minimum Q Q . Region I corresponds to solutions with < m q < 0.91. In this range, 
for each m q there is only one regular D7 embedding, reaching the horizon. In the phase transition 
region (0.913 < m q < 0.926), there are three embeddings for each m q (see the text for details). 
Of the two that reach the horizon one lies on the straight line in region II, the other in region 
III (dashed line). 8 The solution that does not touch the horizon has y m %n = Qo = 0, thus lying 
at the origin IV (we use a dot to describe this kind of solutions). Finally for large enough quark 
masses (m q > 0.93), there is again only one regular embedding for each m q , which does not 
reach the horizon. These solutions all accumulate in the origin IV. 



z (y) V(Q) Q 




Figure 5: Solutions in critical region: There are three condensate values which give different 
regular solutions for the each quark mass. m q = 0.918 in the plot with condensate values 
d = -0.0165, c 2 = -0.0172, c 3 = -0.0265. 



8 There are two small gaps in the plot, one between region II and III and one between III and IV. It is worth 
stressing that quark masses for solutions in these areas all lie in two very narrow ranges around m q ~ 0.925 and 
m q ~ 0.916. At a certain point the fine-tuning required to find the solutions which should fill the gaps exceeds 
the numerical accuracy used to solve the equation of motion. 
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It was known previously that this theory has a first order phase transition at a critical 
temperature of order of the quark mass, and that the chiral quark condensate (ipip) varies 
discontinuously across the transition [201 1221- By computing the effective potential on a 
mixed Coulomb-Higgs branch, we have found that the Higgs VEV is an order parameter for this 
transition. The VEV is zero for < T < T c and non-zero for T > T c . For vanishing quark mass, 
in which case a pure Higgs branch with vanishing Coulomb-branch moduli exists, the effective 
potential suggests that the finite temperature theory is in a (stable) Higgs phase. Although there 
are numerous examples of high temperature symmetry breaking in non-supersymmetric theories, 
our strong coupling result differs from the typical weak coupling behavior of renormalizable 
supersymmetric theories with flat directions, for which the origin of moduli space becomes a 
stable or metastable vacuum at finite temperature. 

Another application of the machinery we have developed here would be to compute the 
effective potential generated on the Higgs branch by other supersymmetry breaking deformations 
of the M = 2 theory, besides the thermal deformations we discussed here. For example, there 
has been considerable interest in deformations that give rise to electrically confining theories. 
In principle, for all such theories one should check that the scalar potential has a minimum at 
the origin of the Higgs branch, before computing physical observables such as meson spectra. 
Moreover, in the future it might be possible with these techniques to reproduce the Affleck- 
Dine-Seiberg runaway scalar potential of M = 1 supersymmetric QCD [SU EH], which arises 
when the number of flavors is less than the number of colors. 
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Abstract 

We use AdS/CFT duality to study the thermodynamics of a strongly coupled M = 2 super- 
symmetric large N c SU (N c ) gauge theory with Nf = 2 fundamental hypermultiplets. At finite 
temperature T and isospin chemical potential /i, a potential on the Higgs branch is generated, 
corresponding to a potential on the moduli space of instantons in the AdS description. For 
fi = 0, there is a known first order phase transition around a critical temperature T c . We find 
that the Higgs VEV is a suitable order parameter for this transition; for T > T c , the theory is 
driven to a non-trivial point on the Higgs branch. For fi ^ and T = 0, the Higgs potential 
is unbounded from below, leading to an instability of the field theory due to Bose-Einstein 
condensation. 
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1 Introduction 



In its original form, AdS/CFT duality [1, 2, 3] relates theories of closed strings in asymptotically 
AdS spaces to large N c gauge theories with matter in the adjoint representation. Fields in 
the fundamental representation may be added by including an open string sector through the 
introduction of branes probing the supergravity background. Much effort has gone into studying 
dualities of this type, motivated largely by the goal of finding a supergravity background dual 
to QCD. 

The first example of AdS / CFT duality for a theory with fundamental representations related 
a conformal M = 2 Sp(N) gauge theory to string theory in AdS§ x S 5 /Z 2 , with D7-branes 
wrapping the Z 2 fixed surface with geometry AdS*, x S 3 [4, 5]. A different approach to flavour 
in AdS/CFT has been considered in [6]. In [7], the duality of [4, 5] was extended to an M = 2 
SU(N C ) theory with Nf massive fundamental hypermultiplets, essentially by removing the Z 2 
orientifold, which was justified by the fact that a probe D7-brane wrapping a contractible S 3 does 
not lead to a tadpole requiring cancellation. Although the field theory is not asymptotically free, 
it has a UV fixed point in the strict N c — > oo limit. Following this, there have been a number of 
papers generalizing the duality to confining theories with fundamental representations [8] -[19], 
including non-supersymmetric examples in which spontaneous chiral symmetry breaking by a 
ijjijj quark condensate occurs [20] -[29]. 

At zero temperature and vanishing quark mass, the M = 2 gauge theory has a non-trivial 
Higgs branch, i.e. a moduli space of vacua on which the scalar components of the fundamental 
hypermultiplets have expectation values. For finite quark mass, the moduli space includes a 
mixed Coulomb-Higgs branch: the fundamental hypermultiplets may have non-zero expectation 
values if the vector multiplet scalars have particular VEV's equal to the quark mass. For both 
the massless and massive case, the moduli space of field theory is described in the AdS picture 
by instanton configurations on the D7-branes [30, 31, 32, 33]. Self dual field strengths are 
solutions of the D7-brane equations of motion due to a conspiracy between the Yang-Mills 
and Wess-Zumino terms in the D7-brane action. In general, this conspiracy is destroyed by 
supersymmetry breaking deformations, which give rise to a potential on the moduli space of 
instantons, corresponding to an effective potential on the Higgs branch. We will compute 
this potential at finite isospin chemical potential and temperature. We focus on the potential 
generated on a slice of the moduli space corresponding to a single instanton centered at the 
origin. This slice is parameterized by the instanton size, which is dual to a particular Higgs 
VEV. 

At finite isospin chemical potential, we find the expected result that the theory is destabilized 
via an effective negative mass squared for the moduli, leading to Bose-Einstein condensation. 
The negative mass squared term in the effective potential is related to the metric on the Higgs 
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branch, which is correctly reproduced by the dynamics of spinning instantons on the D7-branes. 
To stabilize the theory, a positive scalar mass which is larger than the chemical potential would 
need to be introduced. 

Next we consider finite temperature deformations. Some finite temperature properties of 
theories with fundamental representations have been studied previously using AdS/CFT duality 
[20, 22, 34]. As shown in [20, 35], there is a finite temperature first order phase transition as 
the ratio of the temperature to the quark mass is varied in the M = 2 theory of [7] (a similar 
transition was discussed in [22]). In the AdS description, this transition corresponds to a change 
in the topology of the D7-brane embedding in an AdS-Schwarzschild background. 

We explore the finite temperature behavior of this theory in more detail by computing the 
effective potential generated on the Higgs component of the moduli space. We find that the 
Higgs expectation value is an order parameter for the first order phase transition. At < T < T c 
and n — 0, the Higgs VEV is driven to the origin of moduli space. However for T > T c we find 
a surprise; the instanton size is driven towards a non-zero value, suggesting the existence of a 
vacuum in which the theory is higgsed. This is in contrast to the weak coupling behavior of the 
theory, for which the one-loop finite temperature effective potential implies that the origin of 
moduli space is at least metastable. 

The organisation of this paper is as follows. In section 2, we review the AdS description 
[30, 31] of the mixed Coulomb-Higgs branch in the Karch-Katz M = 2 theory [7]. In section 
3, we consider the same theory at finite chemical potential. Section 4 is devolved to the finite 
temperature case. In section 5 we briefly summarize our results and discuss future developments. 

2 SUGRA dual of an J\f = 2 theory with fundamental 
representations 

We consider an M = 2 gauge theory which is dual [7] to string theory in AdS§ x S* 5 with Nf D7- 
branes wrapping a surface which is asymptotically AdS^ x S 3 . The matter content of this gauge 
theory is that of the M = 4 SU(N C ) gauge theory together with Nj massive hypermultiplets in 
the fundamental representation. In Af = 1 superspace, the Lagrangian is 




) 



(i) 



where the superpotential W is 



W = tT(e IJK $i$j$K) + Q l (m + $ 3 )Q l . 



(2) 
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Coordinates 
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D3 
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Table 1: Index conventions 



The superfields Q % and Q iy labeled by the flavor index i = 1---Nf, comprise the Af = 2 
fundamental hypermultiplets. 

This theory is not asymptotically free and, at finite N c , the corresponding string background 
suffers from an uncancelled tadpole. However, as in [7], we focus strictly on the N c — > oo limit 
with fixed Nf. In this case there is a non-trivial UV fixed point for the 't Hooft coupling, 
while the dual AdS string background does not suffer from a tadpole problem since the probe 
D7-branes wrap a contractible S 3 . 

In coordinates which will be convenient for our purposes, the AdS 5 x S 5 background is 

2 p2 ^ 

d, s 2 = ^—dx^dx" + — (dy 2 + y 2 dtt 2 3 + V d^dZ*), (3) 
Rr r z / — ' 

i=l 

e = g s , 

4 

F (5) = dCu) = 4R 4 (V S 5 +* V s5 ), C (4 ) 1 0123 = dx° A dx 1 A dx 2 A rfx 3 , 

it 

R 4 = 4irg s N c a' 2 , 

where is the volume form on S 5 , and r 2 = y 2 + Z % Z % . The 't Hooft coupling in the dual 
gauge theory is A = g 2 N c = g s N c . 

The fundamental hypermultiplets arise from Nf D7-branes embedded in this geometry. Their 
action is 

S D 7 = ~T 7 J d 4 xd 4 y^-detG ab \ PB = - T 7 j d 4 x dy dtt 3 y^l + {Z' i )\ (4) 

where Gpb is the pull-back of the metric, a prime indicates a derivative with respect to y, and 
we have assumed an embedding independent of the coordinates x as well as the coordinates on 
S 3 . There is a manifest minimum of the action when Z[ = 0. We will choose 

Z 1 =m, Z 2 = , (5) 
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other solutions being related by the U(1)r symmetry corresponding to rotations in the plane 
spanned by the Z l . The induced metric is 



r 2 R 2 
ds 2 = —dx 2 / + — (dy 2 + y 2 dQ 2 3 ), 



R 2 



(6) 



with r 2 



ir 



m 



The parameter m corresponds to the mass of the fundamental hypermul- 



tiplets. For m — 0, the geometry is AdS§ x S 3 , while for m ^ 0, the geometry approaches 
AdS 5 x S 3 at large r. The S 3 component of the D7-geometry contracts to zero size at r = m. 
So long as Nf is held fixed in the limit N c —> oo with fixed A = g s N c ^> 1, one can neglect the 
back-reaction of the D7-branes on the bulk geometry. 



2.1 The Higgs branch 

When the theory has massless quarks, the fundamental scalars q % and ^ (denoting bottom 
components of chiral superfields by lowercase letters,) have non-zero expectation values on 
the Higgs branch 1 while the adjoint scalar 03 of the TV = 2 vector multiplet vanishes. For 
non-zero mass, m, and vanishing 3 , the fundamental hypermultiplets are massive and a pure 
Higgs branch does not exist. However there is a mixed Coulomb-Higgs branch when <fr 3 has an 
expectation value such that some components of the hypermultiplets are massless. An example 
of a point on a mixed Coulomb-Higgs branch is given by a diagonal 03 for which all but the last 
k entries are vanishing, 



/o 



\ 



-m 



(7) 



-771/ 



lr rhe scalars 4>\ and fa belonging to the adjoint hypermultiplet may also have non-zero expectation values 
on the Higgs branch. 
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In this case, the F- flatness equations gi(03 + m) — (03 + m)q l = permit fundamental hyper- 
multiplet expectation values in which only the last k entries of q l and q~i are non-zero, 



ft 



(0 



(3 



li 



Pki) 



There are additional F and D-flatness constraints which we have not explicitly written. 

In string theory, nonzero entries in (8) physically correspond to D3-branes which are coinci- 
dent with and dissolved within the D7-branes. Dissolved D3-branes can be viewed as instantons 
in the eight-dimensional world- volume theory on the D7-branes [36], due to the Wess-Zumino 
coupling 

S wz = ^(2W) 2 / C\f B A tr(F A F) . (9) 

There is a one to one map between the moduli space of Yang-Mills instantons and the Higgs 
branch of this M = 2 theory. The ADHM constraints from which instantons are constructed 
[37] are equivalent to the F and D-flatness equations [38, 39] (see also [40] for a review). 

2.2 Supergravity description of the Higgs branch 

Because of the known one-one correspondence between instantons and the Higgs branch, one 
expects that instantons solve the equations of motion of the non-Abelian 2 Dirac-Born-Infeld 
action describing D7-branes embedded in (3) according to (5). The existence such solutions is 
a non-trivial consequence of AdS/CFT duality [30, 31, 41]. 

The effective action describing D7-branes in the AdS background (3) is 



,5' 



T 7 (2wa 



>\2 



— / d 8 £C (4) Atr(FAF) + / d 8 £ ~^-V- det G Tr (F a ^F a(3 ) 
9s J J ^ 



(10) 



where we have not written terms involving fermions and scalars. This action is the sum of a 
Wess-Zumino term, a Yang-Mills term, and an infinite number of corrections at higher orders 
in a' indicated by • • • in (10). The correspondence between instantons and the Higgs branch 
suggests that the equations of motion should be solved by field strengths which are self-dual 
with respect to a flat four-dimensional metric. 



2 The existence of a Higgs branch requires at least two flavors, or two D7-branes. 
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In this paper, we work to leading order only in the large 't Hooft coupling expansion gener- 
ated by AdS / CFT duality, which allows us to neglect the higher order terms in the a' expansion 3 
of the action. Constraints on unknown higher order terms arising from the existence of instanton 
solutions, as well as the exactly known metric on the Higgs branch, were discussed in [30, 41]. 

The induced metric (6) can be written as 



r 2 R 2 4 

ds 2 = —dx^dx 11 + d y md y m > ( u ) 

m=l 

with r 2 = y m y m +m 2 . Field strengths which are self dual with respect to the flat four- dimensional 
metric ds 2 = Ylm=i dy m dy m solve the equations of motion, due to a conspiracy between the 
Wess-Zumino and Yang-Mills term. Inserting the explicit AdS background values (3) for the 
metric and Ramond-Ramond four-form, into the action (10) for D7-branes embedded according 
to (5), with non-trivial field strengths only in the y m directions, gives 



/ 2 f 4 4 r4 f 1 \ 

S Tf(2,7l(y ) I d xd y — I ~Z^mnrsF m nFrs ~\~ F mn F mn I 

J ^gsft V 1 J 

= -T t (4to') 2 / d A xd A y 



4 



(12) 



where F~ n = \{F mn — \e mnrs F rs ). Field strengths F~ n = 0, which are self-dual with respect to 
the flat metric dy m dy m , manifestly solve the equations of motion 4 . These solutions correspond 
to points on the Higgs branch of the dual M = 2 theory. If m ^ this is, strictly speaking, 
a point on the mixed Coulomb-Higgs branch, with expectation values of the form (7), (8). 
We emphasize that in order to neglect the back-reaction due to dissolved D3-branes, we are 
considering a portion of the moduli space for which the instanton number k is fixed in the large 
limit. 

2.3 A slice of the Higgs branch 

For simplicity, we consider the case Nf = 2, which is the minimum value for which a non-trivial 
Higgs branch exists. We will focus on a slice of the Higgs branch (or mixed Coulomb-Higgs 
branch) corresponding to a single instanton centered at the origin, y m = 0. 

In "singular gauge", the SU{2) instanton (using the same conventions as [40]) is given by 

A -0 A - 2< 2 2 ^" m ^ TO (ION 

A 1 " ' A m- y2{y2 + Q2 y W 



3 The dimensionless expansion parameter is a'/R 2 = l/\/A- 

4 Anti-instantons, with F + = 0, correspond to non-supcrsymmetric configurations which do not solve the 
equations of motion. 
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where Q is the instanton size, and 



cr,. 



fzf,l 



2x2J 



cr,. 



cr 



L = (-ir, 1 2X 2) 



(14) 



with r being the three Pauli-matrices. We choose singular gauge, as opposed to the regular 
gauge in which A n = 2a mn y m / (y 2 + Q 2 ), because of the improved asymptotic behavior at large 
y. In the AdS setting, the Higgs branch should correspond to a normalizable deformation of 
the AdS background at the origin of the moduli space. We see that the solution falls to zero 
at large y and the deformation in the interior is controlled by the dimension two parameter Q 2 
which corresponds to the VEV of q l qi. The choice of singular gauge will be of particular use 
when computing the effects of the chemical potential. The singularity of (13) at y m = is not 
problematic for computations of physical (gauge invariant) quantities. 

For m = and zero instanton number, the background geometry preserves the symmetries 
50(2,4) x SU(2) L x SU(2) R x U(1) R x SU(2) f . The 50(2,4) isometry of the AdS 5 factor 
corresponds to the conformal symmetry of the dual gauge theory. The SU(2)l x SU(2)r ~ 
50(4) acts in the obvious way on the coordinates y m . The SU(2) l factor corresponds to a global 
symmetry, while SU(2)r and U(1)r, which acts on (Z 1 ,Z 2 ), correspond to R-symmetries 5 . 
SU(2)f is a gauge transformation on the D7-branes which is constant at the boundary of AdS. 
The action at the boundary of AdS corresponds to the flavor symmetry of the dual gauge theory. 

In the presence of the instanton (13), and for m^O, the symmetries are broken to 50(1, 3) x 
SU (2)l x diag(5[/(2)# x SU(2) f) corresponding to a point on the mixed Coulomb-Higgs branch 

/ \ /o \ /o \ 







\ 







V 



— m 



(15) 



/ 



where qi a are the scalar components of the fundamental hypermultiplets, with flavor index 

2 = 1,2 and SU(2) R index a = 1,2. The scalars <p a belong to the adjoint hypermultiplet, while 

03 is the adjoint scalar belonging to the vector-multiplet. 

In general, one expects that supersymmetry breaking deformations, such as a finite chemical 
potential or temperature, lift the vacuum degeneracy of the Higgs and Coulomb branch. In 
the dual supergravity description, the lifting of the Higgs branch at finite temperature occurs 
because the conspiracy between the Wess-Zumino and Yang-Mills terms, necessary for self- 
dual field strengths to be solutions, no longer occurs. The potential which is generated on the 
Higgs branch can be computed by evaluating the D7-brane action on the space of self-dual field 



5 Thcre is no apparent distinction between SU{2)l and SU(2)r in the geometry of the AdS background. The 
distinction arises from the gauge field couplings to the RR four-form. 



s 



strengths. In the following we will compute the potential V(Q) generated on the slice of the 
Higgs branch dual to the single instanton (13). 

3 Finite chemical potential, zero temperature 

As the simplest example of a potential generated on the Higgs branch, we first consider the case 
of finite chemical potential and zero temperature. For two flavors, there is a U(2) flavor symme- 
try, with Cartan generators 1 and 03, corresponding to baryon number and isospin respectively. 
We will consider a nonzero chemical potential for the isospin 6 . To include the chemical potential 
[42] we allow a spurious gauge field associated with the r 3 component of isospin to acquire a 
VEV, /1, in its A component. This includes generic fermion and scalar Lagrangian terms for 
fields with isospin charge e of the form 

5C = -^rW + A 2 !^ 2 . (16) 

The first term is a source for the fermionic isospin number density. In the path integral, this 
term places the theory at finite density. The second term is an unbounded scalar potential 
which renders the theory unstable, such that Bose-Einstein condensation is expected. We will 
reproduce this run-away behaviour of the scalar potential using the AdS/CFT description. 
Equivalently, we may set A T = by a gauge transformation e lMer , and perform the path integral 
with boundary conditions corresponding to a spinning configuration space. 

AdS/CFT duality relates global symmetries of the boundary theory to gauge symmetries 
in the bulk supergravity. The prescription for turning on a chemical potential in the AdS 
description is to turn on a background (non-normalizable) flat gauge connection A T = \x for 
the associated gauge symmetry. Equivalently, one may consider a spinning AdS background. 
Examples of the AdS description of a chemical potential in various contexts have appeared in 
[43]-[48]. 

The effect of the chemical potential on the Higgs branch can be studied in the AdS description 
by computing the action for rigidly rotating instantons with moduli M 1 , 

A n = eV^mstanton^ ^y-i^ _ ^ 

Equivalently, we may add a background A in the fixed instanton background, 

A = [V _M , A n = A™ tanton . (18) 

On the slice of the Higgs branch corresponding to the single instanton configurations (13) with 
modulus Q, the effective potential at quadratic order in \x can be determined by inserting (18) 

6 Turning on a chemical potential for the baryon number has no apparent effect here since none of the fields 
on the D7 brane world volume are charged under this symmetry. This was pointed out to us by A. Karch. 
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into the D7-brane action J d 4 x V(Q) = —Sm giving 

V(Q) = T r^^- J Atr (\ {y2 + R f ? F~ n F- n + 2F mfl F m ^ (19) 

with y 2 = y m y m . We have split the action into the pieces involving F in the x and y directions, 
indicated by Greek and Roman indices respectively, as well as mixed terms. For the background 
(18), the only non-zero contribution to the potential comes from the mixed term tr F lim F l/m r] tlL ' = 
-tr [A Q , A n } 2 , giving 

2(4W) 2 2 f JA Q 4 m 2(4vrV) 2 2 ^ 2 



V(Q) = -T 7 ^ '-^ / d 4 y " = -T 7 ^ >-^Q 2 . (20) 

9s J y 2 {y 2 + Q 2 ) 2 9s 

Note that the term ^fg~F lim F [1,m in the D7-action is also the term which determines the 
metric (two-derivative effective action) on the Higgs branch. To quadratic order in the chemical 
potential, the result (20) is exact. Higher dimension operators in the DBI action with two Greek 
indices, which have been omitted from (19), must vanish in the instanton background, otherwise 
there would be corrections to the Higgs branch metric in a strong 't Hooft coupling expansion 
[30, 31]. By the same token, the term in the potential on the Higgs branch which is quadratic 
in the chemical potential is given exactly by the F^ m F^ m term. 

In light of (20), the M = 2 theory is unstable at finite chemical potential. Physically, the 
instability is due to Bose-Einstein condensation. This instability will not be cured by taking 
into account the higher derivative terms in the effective action, whose contribution does not 
effect the large Q asymptotics. Potentially, the instability could be removed by the inclusion of 
an appropriately large mass term for the scalars, either as an explicit additional perturbation, 
or generated radiatively from the inclusion of some other relevant perturbation. The finite 
temperature configurations studied below, however, will not be able to cure this instability 
since - as shown below - the relevant deformation, T 4 , has dimension four and can at best 
generate a potential T 8 /Q A . This does not stabilize the theory. 

4 Finite temperature, zero chemical potential 

Before investigating the large 't Hooft coupling behavior of the effective potential on the Higgs 
branch, let us review the situation to one loop in perturbation theory. In a supersymmetric 
theory, the one-loop effective potential is 

d 2 W 



V = Vo + \t 2 J2 



(21) 
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where Vq is the zero temperature effective potential, W is the superpotential, and 0j are the 
scalar components of chiral superfields. For the theory we are considering, this potential drives 
the system towards the origin of moduli space. In fact it is often said that supersymmetric 
theories with flat directions do not have symmetry breaking vacua at finite temperature [49] 
(for some exceptions, see [50, 51]). However, for large 't Hooft coupling, we will find that the 
theory we consider does indeed exhibit symmetry breaking at high temperature. The origin of 
the Higgs branch becomes a maximum of the effective potential for T > T c . 

The gravitational dual of M = 4 gauge theory at large 't Hooft coupling and finite tem- 
perature is obtained by replacing the background (3) with the AdS-Schwarzschild black-hole 
background [53]. The latter belongs to a general class of supergravity solutions which, in a 
choice of coordinates convenient for our purposes, have the form 

4 2 

ds 2 = f(r)(dx 2 + g(r)dr 2 ) + h{r)(%2 dy m dy m + ^dZ i dZ i ), 



e 



m=l i=l 



F (S) = 4tf(y SB + * v s s) = dC {4) , C (4 )|oi23 = s(r) dx° A dx 1 A dx 2 A dx 3 , 

r 2 = y m y m + Z*Z\ (22) 
For the AdS-Schwarzschild solution, we have 

£f , Ar^ + b* , , /4r 4 -6 4 \ 2 , . . R 2 r 4 / b 8 \ 

/M = • s(r) = \4r*T¥ ) • hir) = T 5 "' = W\ 16? j (23) 

The coordinates x are the spatial coordinates of the dual gauge theory and r is the Euclidean 
time direction, which is compactified on a circle of radius b~ l , corresponding to the inverse 
temperature. Note that the temperature T ~ b only enters to the fourth power. 

The D7 embedding in this background to describe our M = 2 gauge theory with quarks is 

Z 2 = 0, Z x = z{y), (24) 

where y 2 = y m y m . The induced metric G| p b takes the form 

ds 2 m = f(r) (dx 2 + g{r)dr 2 ) + h(r) (^dy m dy m + (^^j (y m dy m ) 2 J , (25) 

with r 2 = y 2 + z 2 (y). To lowest order in a', the D7-brane action is 

-TV y ^ e-^v^detGlpt = -2vr 2 T 7 J d A x dy y 3 0(r)/(r) 2 MOVl + ^'(y) 2 . ( 2 6) 

and the specific embeddings z(y) can be found by solving the Euler-Lagrange equation with 
boundary conditions z(oo) = m, z'(oo) = 0. 
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(27) 



At large y the geometry returns to AdS where the asymptotic embedding solution (24) takes 
the form z(y) ~ m + c/y 2 . Given m, the parameter c, corresponding to a tpip condensate, 
is determined by requiring the solution to be smooth and normalizable. These solutions were 
studied in detail in [20]. 

The potential generated on the Higgs branch is obtained by calculating a contribution to the 
D7-brane action involving the gauge field F a b at second order in a'. Specifically, one evaluates 
the action on the space of field strengths which are self-dual 7 with respect to the induced metric 
in the directions transverse to r, x; 

V = T7(2 ^ )2 (j J d A yC^e mnrs TiF mn F rs - l - J d 4 y V^drtG Tr F mn F mn ^ , 

where F mn is self-dual with respect to the metric 

ds\ = Kr) ((1 + z\yf)dy 2 + y 2 dn 2 3 ) . (28) 

This metric is conformally flat. With new coordinates y(y) such that ds 2 = a(y)(dy 2 + y 2 <if2§), 
the instanton configurations (self-dual field strengths) take the usual form. 

In the b = (zero temperature) limit, the analysis of section (2.2) holds and the potential 
(27) becomes flat. At finite temperature, the conspiracy between the Wess-Zumino term and the 
Yang-Mills terms fails, giving a potential on the moduli space of instantons. We now compute 
the form of the potential for the slice of moduli space corresponding to an instanton centered 
at the origin. 

4.1 Computation of the Higgs Effective Potential 

We wish to evaluate the D7-brane action on the space of fields strengths which are self-dual 
with respect to the transverse part of the metric 

dy m dy m + f^M] {y m dy m " 2 



y 

This metric is conformally flat. There exists a coordinate transformation 

y m = J(y)y m (29) 

such that the induced metric has the form 

dsl, = f(r) (dx 2 + g(r)dr 2 ) + -^dy m dy m . (30) 



7 Thcre are couplings between world-volume scalars and field strengths at higher orders in a' which could 
alter the embedding. However we only consider the leading term in a large 't Hooft coupling expansion for which 
these couplings can be neglected. 
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The function J(y) is obtained by solving 



J'Y 2 J' 1 ,. , 2 

j) +---- z (y) = 



(31) 



y J y 2 

subject to the boundary condition J(y — »• oo) = 1. In the coordinates y m , the single instanton 
centered at the origin has the usual form, 

— ~2(~2 i /02^ ' y ' 

y 2 (y 2 + Q 2 ) 

The embedding function z(y) is determined by shooting techniques as in [20]. The potential on 
the Higgs branch is then given by Sd7 = — f d 4 xV(Q), or 

(2W) 2 



V(Q) = T 7 
= T 7 
= T 7 



2 

^9s 



d 4 y (s{r) - (p{r) f{rf ^/gjr) 



96 g 4 
(y 2 + Q 2 ) 



(33) 



Vol(S 3 ) / 2/ 3 dyfs(r)-0(r)/(r) 2 v^(7 



96 Q 4 J{y) 4 

(J(y) 2 y 2 + Q 2 ) 4 



Ay) 
J(y) 



For the AdS-Schwarzschild background (22,23) we obtain 

6(2vrV) 2 f, 6V Q*J(y) 4 



V(Q) = T 7 - 



g s R 4 



dy 



(J(y)2y2 + z{y) 2 ) 2 { J {y) 2y2 + Q 



2\4 



Ay) 
Ay) 



(34) 



4.2 Large VEV Asymptotics 



Without resorting to numerics, some qualitative features of V(Q) at large Q (or large instanton 
size) can be determined from the large y behavior of the induced metric and RR-form given in 
(22,23), together with the asymptotic behaviour of the embedding, 



z(y) = m + — + 



2T 



(35) 



At large y we can neglect the condensate and the leading terms in the potential will depend 
on m and on the deformation parameter of the background b. The D7-brane action for a large 
instanton, which has support at large y is 



J d 4 xV(Q) = -S D7 



2„,A2 



T 7 



6(2vr 2 a 
9sR 4 

6(2ttV) 2 
9sR 4 



d 4 x / dy y 3 



Q 4 



(y 2 + Q 2 ) 4 



b 8 2b 8 m 2 



y> 



Q 4 



- log 2 - - 

2 & 3 



2m 2 



<> 



y 



67 
48 



+ 



2 log 2 



(36) 
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This asymptotic approximation is valid when y > Q but gives IR divergent integrals. We have 
evaluated the integrals down to y/Q = 1. In the next section we evaluate the full potential 
including the correct, convergent, IR behaviour. 

The asymptotic potential terms are simply what one would expect on dimensional grounds. 
The potential must vanish as b — > and the geometry becomes AdS. Since b enters as b A , the 
first Q dependent term must take the form b 8 /Q 4 . 

4.3 Numerical computation of V(Q) 

To compute V(Q) in general requires knowledge of the embedding function z(y), which has 
been computed by a numerical shooting technique in [20]. Imposing boundary conditions cor- 
responding to the large y behaviour (35), and requiring smooth behaviour in the interior, such 
that an RG flow interpretation is possible, leads to a dependence cim, b) of the chiral quark 
condensate ('tpip) on the quark mass and on the temperature. Depending on the ratio m/b, 
there are two types of solutions, which differ by the topology of the D7-branes. At large r 
(or y) the geometry of the D7-branes is AdS 5 x S 3 and the topology of the r — > oo boundary 
is S 1 x R 3 x S 3 . For sufficiently large m/b, the S 3 component of the D7-geometry contracts 
to zero size at finite r > b. In this case the D7-brane "ends" before reaching the horizon at 
r = b. However, for sufficiently small m/b, the D7-brane ends at the horizon, at which point the 
thermal S 1 contracts to zero size. Both these types of solutions are plotted in figure 1. There is 
a first order phase transition at the critical value of m/b ~ 0.92 where the two types of solution 
meet [20, 35]. The (ipip) condensate is non-zero on both sides of this transition, although there 
is a discontinuous jump in its value. 



z(y) 




m=0.92 
m=0.91 



Figure 1: Brane embeddings in the AdS Schwarzschild background for different values of the 
quark mass. In the plot we have set b = 1. 
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Using the machinery discussed above, we can study the behaviour of the scalar VEV 
parametrizing the Higgs branch on either side of the transition. We evaluate the potential 
(27,34) numerically, with a single instanton centred at the origin for different values of the ratio 
m/b. The results are plotted in figure 2. 




V(Q) V(Q) 



Figure 2: Potential V(Q) as a function of the instanton size / Higgs VEV Q for various values 
of the quark mass m (we set 6 = 1 here); in fig. a) we sample the whole spectrum from m = 
to m — > oo, in fig. b) we vary m close to the phase transition region. The potential for m — » oo 
is flat and coincides with the horizontal axis. 



Figure 2 shows that there are two phases for the scalar Higgs VEV matching the two phases 
(topologies) of the D7-brane embedding. These phases are divided by the first order phase 
transition at a critical value of m/b, where the parameter b of the AdS/Schwarzchild background 
is proportional to the temperature. In figure 2 we have chosen units such that 6=1. For 
m/b — > oo, the D7 branes probe the AdS-like part of the geometry, such that the potential is 
flat (dotted line on the left hand side of figure 2). For all other values of m/b, the potential 
approaches this constant asymptotic value from below for large Q. When the value of m/b is 
decreased, starting at m/b — > oo, a minimum of the potential forms at Q = 0, with increasing 
depth. At the critical value m c /b ~ 0.92, a new solution for the brane embedding appears, for 
which the potential has a minimum at Q ^ 0. For m/b — > 0, V has a maximum at Q = and 
a minimum of decreasing depth for Q ^ 0. Atm/6 = (solid line in the plot) there is still a 
shallow minimum of V at Q ~ 0.856. Thus, the Higgs VEV is an order parameter for the first 
order phase transition. We emphasize that the symmetry breaking phase occurs at small m/b, 
or high temperature, contrary to the usual expectation. We plot the value of the Higgs VEV as 
a function of m/b in Figure 3. 
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Figure 3: Position of the minimum of the potential Qo versus the bare quark mass m q , zoom of 
the critical region. 

Note that the expectation value Qo is a multivalued function of m/b in a region near the first 
order critical point. There are several regular embeddings in this region, with the true vacuum 
corresponding to the solution with lower free energy. 

Thus we find that for T > T c , the theory is driven to a Higgs phase or, more precisely, a point 
on the mixed Coulomb-Higgs branch. We have not determined whether or not this vacuum is 
merely metastable, which would be the case if the origin of the Coulomb branch has lower free 
energy. 

Since the less supersymmetric phase occurs for T > T c , the phase transition discussed here is 
quite different from the deconfinement transition. Nevertheless it is interesting to note that our 
phase transition is first order, as is found in lattice simulations of the deconfinement transition 
in large N SU(N) gauge theories [52]. 

4.4 Geometric aspects of the transition 

It is worth stressing the link between the transition of the Higgs VEV between the two phases 
and the change in the embedding topology of the D7 brane. From figures 1 and 2, we see that 
the Higgs potential only has a non-zero minimum when the D7-branes reach the horizon. There 
is a repulsive effect in the vicinity of the horizon which causes the instanton to expand. In 
figure 4 we plot the value of the Higgs VEV versus the minimum value of the coordinate y to 
which the D7 branes extend. When the D7-branes do not reach the horizon, the S 3 component 
of their geometry contracts to zero size at the endpoint, which occurs at y — 0. On the other 
hand, when the D7-branes reach the horizon, the S* 3 component of their geometry remains finite 
size and y ^ at the endpoint. The dependence of the VEV on the minimum value of y is 
almost a straight line with slope roughly equal to 1.25. The solutions that reach the horizon are 
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distributed along this line, while the ones which end before reaching the horizon accumulate at 
the origin. 

The phase transition region (m/b ~ 0.91) has some additional interesting structure. In this 
region there are three regular solutions for the D7 embedding for each value of m q (see [35, 22]). 
Some of them end on the horizon and some do not. The three flows are shown for a particular 
case in fig. 5a. These different embeddings give different predictions for the Higgs VEV. We also 
show the potential for Q on each of these embeddings in figure 5b and the minimum vs m can be 
found in figure 3b. There are three values of Q supplied in this transition region. Clearly there 
will be a first order transition where the embedding jumps from solutions ending on and off the 
horizon as m is changed, corresponding to a discontinuous transition in the value of the scalar 
VEV. 

5 Conclusions 

We have investigated the thermodynamics of a large N c strongly coupled M = 2 gauge theory, 
using AdS/CFT duality to compute the potential which is generated on the Higgs branch at 
finite temperature and isospin chemical potential. In the AdS description, this potential is a 
potential on the moduli space of instantons. 

We have shown non-perturbatively that the chemical potential destabilizes the Higgs poten- 
tial, giving rise to Bose-Einstein condensation. The 0(fi 2 ) term in the potential can be computed 
exactly, as a consequence of the non-renormalization of the metric on the Higgs branch, which 
is correctly reproduced by instanton dynamics in the AdS description. The instability due to 
the chemical potential is not cured at finite temperature. 

It would be very interesting to introduce a chemical potential for Baryon number in the AdS 
setting. Unlike the chemical potential for isospin, this is known to be very difficult to study 
on the lattice due to a complex fermionic determinant. In the AdS setting, it is also easier to 
introduce a chemical potential for the isospin, as we have done here, although the problems 
arising for baryon number are probably much more tractable than those which arise on the 
lattice. 
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Figure 4: Position of the endpoint y m in of the embedded brane on the horizon versus the potential 
minimum Qo- Region I corresponds to solutions with < m q < 0.91. In this range, for each m q 
there is only one regular D7 embedding, reaching the horizon. In the phase transition region 
(0.913 < m q < 0.926), there are three embeddings for each m q (see the text for details). Of 
the two that reach the horizon one lies on the straight line in region II, the other in region III 
(dashed line). 8 The solution that does not touch the horizon has y min = Q = 0, thus lying at 
the origin IV (we use a dot to describe this kind of solutions). Finally for large enough quark 
masses (m q > 0.93), there is again only one regular embedding for each m q , which does not 
reach the horizon. These solutions all accumulate in the origin IV. 
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y 

Figure 5: Solutions in critical region: There are three condensate values which give different 
regular solutions for the each quark mass. m q = 0.918 in the plot with condensate values 
d = -0.0165, c 2 = -0.0172, c 3 = -0.0265. 



8 Thcre are two small gaps in the plot, one between region II and III and one between III and IV. It is worth 
stressing that quark masses for solutions in these areas all lie in two very narrow ranges around m q ~ 0.925 and 
m q ~ 0.916. At a certain point the fine-tuning required to find the solutions which should fill the gaps exceeds 
the numerical accuracy used to solve the equation of motion. 



18 



It was known previously that this theory has a first order phase transition at a critical 
temperature of order of the quark mass, and that the chiral quark condensate (ipip) varies 
discontinuously across the transition [20, 35, 22]. By computing the effective potential on a 
mixed Coulomb-Higgs branch, we have found that the Higgs VEV is an order parameter for this 
transition. The VEV is zero for < T < T c and non-zero for T > T c . For vanishing quark mass, 
in which case a pure Higgs branch with vanishing Coulomb-branch moduli exists, the effective 
potential suggests that the finite temperature theory is in a (stable) Higgs phase. Although there 
are numerous examples of high temperature symmetry breaking in non-supersymmetric theories, 
our strong coupling result differs from the typical weak coupling behavior of renormalizable 
supersymmetric theories with flat directions, for which the origin of moduli space becomes a 
stable or metastable vacuum at finite temperature. 

Another application of the machinery we have developed here would be to compute the 
effective potential generated on the Higgs branch by other supersymmetry breaking deformations 
of the M = 2 theory, besides the thermal deformations we discussed here. For example, there 
has been considerable interest in deformations that give rise to electrically confining theories. 
In principle, for all such theories one should check that the scalar potential has a minimum at 
the origin of the Higgs branch, before computing physical observables such as meson spectra. 
Moreover, in the future it might be possible with these techniques to reproduce the Affleck- 
Dine-Seiberg runaway scalar potential of M = 1 supersymmetric QCD [54, 55], which arises 
when the number of flavors is less than the number of colors. 
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